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The oscillation periods of the interlayer exchange coupling are investigated when two magnetic 
layers are separated by a metallic superlattice of two distinct non-magnetic materials. In spite of 
the conventional behaviour of the coupling as a function of the spacer thickness, new periods arise 
when the coupling is looked upon as a function of the number of cells of the superlattice. The 
new periodicity results from the deformation of the corresponding Fermi surface, which is explicitly 
related to a few controllable parameters, allowing the oscillation periods to be tuned. 
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I. INTRODUCTION 

Oscillatory interlayer exchange coupling in metallic magnetic multilayers causes the magnetizations of neighbouring 
magnetic layers separated by non-magnetic spacers to be spontaneously aligned ferromagnetically or antiferromag- 
netically depending on the thickness of the spacer material. This phenomenon has been intensively studied both 
experimentally and theoretically over the last few yearsau. 

The relation between the periods, by far the most investigated feature of the oscillations, and the electronic structure 
of multilayered systems is a fundamental aspect of the mechanism responsible for the origin of this phenomenon. In 
fact, the periods were originally shown to be dependent on the shape of the spacer Fermi surface, YjteJ e two distinct 
criteria were proposed for determining the periods with which the interlayer coupling oscillatcsQBO. Within the 
quantum well theory the wave vectors yielding the effective oscillation periods correspond to extremal radii of the 
Fermi surface in the direction perpendicular to the layers, or in other words to half the caliper measurements! 
On the other hand, in the RKKY theory adapted to the layered geometry, the periods are associatedpivith wave 
vectors perpendicular to the layers linking two points of the Fermi surface with antiparallel velocities™!]. Despite 
the distinction between those geometrical criteria, the periods predicted by both theories, at least in their original 
formulations, agree in most cases. When they do not coincide, such as, for example, when the lattice lacks reflection 
symmetry about the layer planes, suitable transformations of the Fermi, surface can be performed, which reestablish 
the correspondence between the quantum well and the RKKY periodso. Recently howe^eiv-fundarnental oscillation 
periods not directly associated with the shape of the spacer Fermi surface were obtainedB&EI Those periods, which 
correspond to higher order contributions to the exchange coupling, are no longer determined by the above geometrical 
criteria but by a set of more general conditions on the electronic structure of the system, although they are still 
indirectly dependent on the shape of the Fermi surface. These conditions provide a systematic way of determining 
the oscillation periods involved in the interlayer coupling. 

In order to test the applicability of those general conditions, it is worth investigating them in situations where the 
spacer Fermi surface, clearly a relevant factor in this matter, is variable. One approach is to use an alloy spacer of 
variable composition although the Fermi surface is no longer sharply defined. Also, it is well know that the scattering 
of electrons by random impurities strongly damps the oscillations. In fact, the amplitude of the coupling is found 
to decay exponentially with the concentration of impuritiesliiH13. An alternative to an alloy is a metallic superlattice 
whose Fermi surface depends on the proportions and layer thicknesses of the constituent materials but where the 
oscillations are not damped by incoherent electronic scattering. Moreover, nature provides us with only a limited 
number of non-magnetic metals. However, by combining distinct materials in superlattices the number of possible 
non-magnetic superlattices is infinite. 

The purpose of this communication is to investigate the oscillation periods of the interlayer exchange coupling 
when the spacer is a non-magnetic superlattice. It is shown that the alteration of the oscillation periods reflects the 
deformation of the corresponding Fermi surface only when the coupling is investigated as a function of the number 
of superlattice cells separating the magnetic layers. In the simple model treated here, the oscillation periods of the 
coupling are written in terms of controllable parameters which are related to the composition of the superlattices, 
and which, accordingly, allow the periods to be tuned. It is hoped that such a simple model reflects what occurs in 
realistic systems, and this may be confirmed by future experiments. 
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The sequence adopted in this work is described as follows. In order to make the work self-contained, a brief 
summary of the general method for determining the oscillation periods in a multi-orbital tight-binding modelj is 
presented in section |l[ Secondly, the method is shown to be applicable when the spacer is a superlattice, using the 
clear correspondence between orbital and intraccll indices. Subsequently, the oscillation periods are obtained following 
the analysis of the electronic structure of the system and are explicitly written in terms of parameters associated with 
the superlattice composition. 



II. OSCILLATION PERIODS 



A brief summaryi-pf the general method for determining the oscillation periods of the interlayer coupling is presented 
here. See referenced for a comprehensive discussion of the method. 

For the sake of simplicity, we consider two parallel magnetic planes, labelled and n, embedded in an infinite 
non-magnetic metal. As far as the coupling as a function of the spacer thickness is concerned, theunumber of magnetic 
planes influences only the phase and amplitude of the oscillations and does not affect the periodsE3. Since the periods 
are the main concern here, the restriction to two single magnetic planes does not pose any limitations on the results 
obtained. 

The interlayer coupling J is given by the difference in total energy between the antiferromagnetic and ferromagnetic 
configurations or, in other words, by the energy necessary to flip the magnetization in one of the planes. Based on the 
formalism introduced by d'Albuquerque e Castro et al, the coupling J is written in terms of the one-electron Green 
functions of the multilayered system by the following expression 

J = J2 I rf«^(q||,w) . (!) 
ill J 

where JF(qy , a;) = f(u>) ImTr In |l + 4G^ (q||, u>) V ex Gg n (qy, oj) V ex \- In the equations above G^ (q\\,uj) is the 

propagator between planes and n for electrons with spin a in the ferromagnetic configuration of the system, f(ui) is 
the Fermi function, V ex is a diagonal matrix in orbital indices representing the strength of the local exchange potential 
in the ferromagnetic layers, and the sum over qy is restricted to the two-dimensional Brillouin zone. Bearing in mind 
that the exchange splitting of the magnetic material is localized on the two planes and n, the propagators G can be 
rewritten as a function of the bulk spacer propagators g n g and g$ n , i.e., 

Gn (qi . w ) = Qa (Pno(q|| , w), 5on(q|| , w)) . (2) 

Moreover, it can be shown that within the single-band model the dependence of G a upon those propagators is such that 
G^Q = G" (g n o gon)- It is clear that the spacer thickness dependence of the interlayer coupling is entirely determined 
by the off-diagonal propagators g n o and go n . The investigation of the oscillation periods thus requires an analysis of 
g n o as a function of the interplane distance. 

In the layered geometry, where qii represents the in-plane wave-vector, the expression for a general matrix element 
9em representing the propagation from an orbital [i at plane £ into an orbital v at plane m of an electron at the Fermi 
level Ep is given by 

r >_,< i >^ r KMi'-M'-" 1 - ( '-' ,d 



*«ii.ft> - ( 5 ) L l, * X-a( q ,^> + i°+ • (3) 

where the integration variable q± is the wave vector perpendicular to the layers, d is the interplane distance, E s (qii , q±) 
describes the bulk spacer band structure, s is the band index, and a sfl (q±) = (q±s\q±fi) is the projection of the 
eigenstate s into the orbital /i for a given q± . By changing the integration contour from the real axis to the boundaries 
of a semi- infinite rectangle in the complex plane whose base lies on the real axis between —ir/d and ir/d, the evaluation 
of the integral above is simplified to the determination of the poles of the integrand inside the appropriate contour 
and their respective residues. The rectangle is in the upper-half plane for £ < m and in the lower-half plane otherwise. 
Clearly, the poles are given by the values of q± satisfying the condition Ep + — E s (q\\ , q±) = 0, where Ep + = Ep+i + . 

By adding the residues associated with the contributory poles, labelled q?, an analytical expression for gp m is 
obtained, i.e., 



gtrn(<l\\ 



^Aye"^-" 1 ^ , (4) 
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where the matrix elements of the matrix A S j 



d£ s (q||,g_L) 



(5) 



The coefficients A^j depend neither on I nor on m and the only dependence on the distance between the planes is 
inside the argument of the exponentials. It should be stressed that this simple method for obtaining the one-electron 
propagators is exclusive to the one-dimensional case, and therefore for fixed values of . The exponentials in Eq. (Q) 
are independent of the orbital indices, which means that all matrix elements oscillate with the same periods. Note 
that in calculating Eq.(|^) it was assumed that the coefficients a sfJl (q±) have neither singularities nor branch points 
inside the integration contour. In this way, the coefficients are evaluated just at the poles q s j. In general, branch cuts 
associated with the coefficients a sfJ- (q±) may exist, but because they may be chosen not to intersect the real axis, 
the factor e -*u(*-m)d in the integranc j f Eq .(g) 

ensures that the corresponding contribution is negligible for large 
\i — m\. Therefore, the expression above is correct in the asymptotic limit of large interplane distanced. 

It is clear from Eq. (||) that gf m is an oscillatory function of the distance between the planes and it oscillates with 
different periods, each one associated with its respective wave vector g|. For this reason it is interesting to look at the 
physical significance of the poles. The poles may be real or complex, but we note that only the real poles contribute 
significantly to the coupling. This is because complex poles, when inserted in Eq.(^), damp the oscillations and are 
important only in cases of very thin spacers. Real poles, which are obtained from the electronic structure of the bulk 
spacer, are just the perpendicular coordinates of the Fermi surface for a fixed value of qii . In other words, they are 
the perpendicular components of the wave vectors with which electrons of energy Ep propagate across the spacer. 

Having investigated gi m as a function of the interplane distance, it is worth stressing that the functions T must 
oscillate with the-,sanae periods. When expanded in a Fourier series, taking into account the possibility of quasi- 
periodicity of <7 n j3oEj, the function T is written as 

F= £ C (m} ,... ira . i ...,e iwi ^ (m ^ ) > (6) 
(mj,...,mj,...) 

where C( m i m = ) are the Fourier coefficients determining the amplitude of the oscillations and depend on the exact 
form of the function T . The indices [m\, m|, ...) are integers and there are as many as the number of poles q?. Note 
that the periods are given by ^( m \,...,m s ,...) = 2tt/ (J2 S j m j1j) anc ^ do n0 * depend on the exact form of T, since they are 
obtained from the bulk spacer band structure. Because the indices run over the set of integers, there is always an 
infinite number of periods defined by Eq. (^|) , even for a small number of poles q| . In spite of the somewhat congested 
notation of the general expression above, in practice there are only a few indices involved due to the restricted number 
of real solutions crossing the Fermi level. In fact, the number of solutions is related to the number of sheets of the 
Fermi surface. 

When Eq.(|(]) is inserted into Eq. ([[]), the two-dimensional integrali-can be evaluated through the stationary phase 
approximation, valid in the asymptotic limit of large spacer thicknessu'E-3. In doing so, the periods of the coupling are 
identified by selecting the periods of T associated with wave vectors satisfying the following condition, 

2["*jV||^(q„,^)] = 0, (7) 

SO 

where Vn is the two-dimensional gradient in q||-space. The equation above gives the necessary conditions for con- 
structive interference between electrons across the spacer. The otherwise destructive interference does not contribute 
to the coupling. Eq.(0) simply selects the points across the two-dimensional Brillouin zone which yield the effective 
periods with which the coupling oscillates. 

As pointed out in referenced, although the wave vectors are directly related to the Fermi surface, the general 
wave vectors Yl s j( m j Qj) are n °t- This is a result of the interference between incommensurate wave vectors giving 
rise to a new set of periods not directly obtained from the structure of the Fermi surface. Consequently, geometrical 
criteria cannot be used in general to select the oscillation periods. However, for Fermi surfaces with a single sheet, 
where the off-diagonal propagators oscillate as a function of the interplane distanegjpth only one period, the selection 
rules given by Eq.(Q) coincide with the geometrical criteria of the RKKY theorytm. The only difference in this case 
is that Eq.(|^) determines the respective harmonics, in addition to the fundamental RKKY periods. 

In summary, the determination of the oscillation periods results from the analysis of the bulk spacer off-diagonal 
propagators where the relevant real wave vectors, under the stationary phase condition, are obtained from the band 
structure. When dealing with multi-sheet Fermi surfaces, the possibility of interferences between incommensurate wave 
vectors may introduce new features in the oscillations, which are not present in the case of single-period oscillations 
of the propagators. 
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III. SUPERLATTICED SPACERS 



Having summarized the method for determining the oscillation periods of the coupling, it is worth stressing that no 
particular form for the spacer material was assumed. The only information required is the bulk band structure which 
makes possible the investigation of the coupling across arbitrary spacers. Spacers consisting of superlattices of distinct 
non-magnetic materials will be investigated here. More specifically, the material into which the magnetic planes are 
embedded is composed of two non-magnetic metals A and B with respective thicknesses Na an d Nb periodically 
alternated, as schematically displayed in figure |l|. In order to make a distinction, these structures will be hereafter 
named superlatticed spacers, in contrast to the conventional single-layered spacers. The interlayer coupling across 
superlatticed spacers can be looked upon as a function of different thicknesses, namely Na, Nb, N and N c , where N c 
represents the number of cells separating the magnetic planes. Each cell is formed by a set of Na planes of metal A 
and Nb planes of metal B. As far as the oscillation periods are concerned, the coupling across superlatticed spacers 
as a function of Na or Nb does not present new features when compared to the coupling across conventional single- 
layered spacers as a function of Na or Nb, respectively. In other words, if a single-layered spacer of material A yields 
an interlayer coupling with period Xa and likewise As is the period for the coupling across a single-layered spacer of 
material B, the coupling across a superlatticed spacer of materials A and B will oscillate with period A a when seen 
as a function of Na and with period As when seen as a function of Nb- However, when looked upon as a function 
of N c , for fixed values of Na and Nb, the coupling oscillates with periods different from A^ and Xb- Although the 
superlattice Fermi surface is different from the individual Fermi surfaces of materials A and B separately, it is only 
when investigated as a function of N c that the deformation of the Fermi surface is translated into a change in the 
periods. 

Following the steps of the previous section, the bulk spacer band structure must be calculated. For the sake 
of simplicity, the single-band tight-binding model is considered assuming a simple cubic lattice structure for the 
superlattice grown along the (100) direction. For the two given materials A and B, the in-plane energies are labelled 
tA and tB, and the interplane hoppings are tA and ts- The hopping between planes A and B is represented by 
tAB- We recall that the tight-binding parameters above depend on the in-plane wave vector qii through e^(qii) = 
e° + 2t® [cos(q x d) + cos(q y d)} and ^(qy) = tg, where [3 is cither A or B, ejg is the on-site energy and tp is the hopping 
between nearest-neighbour sites in material (3. Note that q z is perpendicular to the layers whereas q x and q y are the 
coordinates of qn . 

Because each set of N planes is periodically repeated, where N = Na + Nb, the system can be represented by cells 
of N planes equally spaced. Consequently, the tight-binding Hamiltonian of the infinite multilayered system, when 
written in the appropriate basis, becomes a N x N matrix given by 

[H A (N A )] [Has] 



U ~\ [H B a] [Hb(Nb)} ) ' (8) 

where [Ha(Na)] and \Hb{Nb)\ are square matrices of order Na and Nb, respectively. The sub-matrices [Ha(Na)] 
and [Hb(Nb)} are explicitly given by 
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> N/3 . (9) 



Note that \Hp(Np)] is a tri-diagonal matrix equivalent to the single-band tight-binding Hamiltonian of a slab consisting 
of Np planes. The sub-matrices off the diagonal in Eq.(||) are Na X Nb and Nb X Na matrices such that 

ooo e - iq±Nd 

[H A b] =^b( j } N A (10) 
10 



and [Wba] = ("Hab) ■ In the matrix directly above q± is the wave vector perpendicular to the layers and d is the 

interplane distance, which is assumed to be the same for both materials. The basis producing such matrices results 
from a Fourier transform of the localized- atomic-orbital basis along the direction perpendicular to the layers. Since 
the cells of N planes are equally spaced by Nd, the associated wave vector q± is in the range [jni, wi], that is, the 
corresponding Brillouin zone along this direction is reduced by a factor N . The tight-binding Hamiltonian above 
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is equivalent to the TV-band Hamiltonian of a one-dimensional system. In other words, the N intracell planes may 
be considered as orbitals of a linear chain whose nearest neighbour elements are Nd apart. Having established this 
correspondence, hereafter intracell and orbital indices will be mentioned with no distinction between them. 

Following section |l[ the oscillation periods result from the analysis of the bulk spacer off-diagonal propagators as 
a function of the distance between the planes. Note however that in the case treated here the effective interplane 
distance must be the distance between the cells, i.e., D = Nd. Hence, the off-diagonal propagators as a function of the 
number of cells N c are the ones to be investigated. These propagators depend on the band structure of the system, 
as clearly shown in Eq.(^). The original tight-binding Hamiltonian is partially diagonalized by both in-plane and 
perpendicular Fourier transforms, the remaining problem being the determination of the eigenvalues of the N x N 
matrix in Eq.(||). Therefore, the information on the oscillation periods of the coupling in a superlatticed system 
depends on the diagonalization of this matrix, which will be presented in the next section. 

To maintain the correspondence between the intracell planes of the superlattice and orbital indices of the linear 
chain, one has to assume that in the case of superlatticed spacers the magnetic planes embedded in the non-magnetic 
system must be formed by two separate magnetic cells. However, for the sake of simplicity and with no loss of 
generality, we consider only two magnetic planes. If we label the magnetic planes with the intracell indices a in the 
cell and a' in the cell N c , the relevant Green functions to be calculated are the off-diagonal propagators ggjv (l|| > w )> 
i.e., 

In addition to the fact that a and a' are intracell indices, the only difference between the equation above and Eq.(^) 
is the interplane distance D which modifies the integral limits. Similarly to section O, the integral above is evaluated 
through the residues method, where the determination of the poles follow the diagonalization of the Hamiltonian. 



A. Electronic structure 



Making use of the particular symmetry of the matrices displayed above, the determinant of TL in Eq.(H) can be 
written in terms of the sizes of its sub-matrices Na and Nb yielding a simple analytical expression for the oscillation 
periods of the one-electron propagators. The energy eigenvalues u are solutions of the equation det \H — lu1\ = 0, 
where 1 is the identity matrix. In this case the sub-matrices along the diagonal become 



[Hp(Np)-.i] = *g * « , (12) 
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whereas the block matrices off the diagonal are still given by [Has] and [Hba]- 

The determinant of a N x N square matrix is a sum of N\ terms. Those terms are made of products of N matrix 
elements, following all their possible permutations. Because the sub-matrices on the diagonal have only two distinct 
non-zero elements, the evaluation of their determinants by counting the number of combinations allowed by the 
possible permutations is not particularly difficult and can be expressed in term of the sizes Np. The solution of this 
combinatorial problem is given byliS 

N f>/ 2 / X 

det\Hp(Np)-ui\ = £ (-iy ( N e7 J ) (e p ~u) N ^ (tp)* , (13) 

where j is an integer, the upper limit of the sum is the integer part of Np/2 and 

fNp-j\ (Np-j)\ 



V 3 J (Np-2j)lj\ (M) 

is the binomial coefficient. Note that the expression above is valid for symmetric tri-diagonal matrices with similar 
elements along each of the diagonals. Bearing in mind that Ttp(Np) is equivalent to the slab Hamiltonian, the values 
of uj satisfying det \Hp(Np) — uj1\ =0 through Eq.(|l3|) are just the associated eigenvalues. 
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The Hamiltonian TL in Eq.(||) is almost tri-diagonal. The matrix TL is not tri-diagonal due to two of its elements, 
namely TLi^n and TLn,i, which arc non-zero. Therefore, an expression for the determinant of TL is not expected to be 
very different from the one in Eq. (|l3|) . In fact, by considering the only two non-zero elements of the submatrices in 
Eq.(fTo|) and recounting the number of possible permutations, the determinant of the matrix TL can be expressed in 
terms of the determinants of the matrices TLa(Na) and TL B (N B ). In the similar case of TL — u>l, the determinant is 
given by 

det \H - wl| = dot \Ha(N a ) - ul\ x det \H B (N B ) - ul\ 

- 2 (t AB f det \H A (N A - 1) - u;i| x det \H B (N B - 1) - wi| 
+ (t AB f det \H A (N A - 2) - wi| x det \H B (N B - 2) - wi| 

+ {t AB ) !> (t A f Ji - 1 (tB) NB - 1 coa(q ± Nd) . (15) 



Note that the equation above is written in terms of lower-order determinants of Eq. ( |13| ) and is an analytical expression 
for det \TL — u>\\ in terms of its sizes Na and N B . 

It is instructive to point out that when A = B, that is, when a single-layered spacer is considered, the usual 



eigenvalues u>{q±) = e + 2tcos(q±d) of an infinite linear chain are obtained from Eq.(15) through the condition 
det|W-wl| = 0. 

The eigenvalues ui of the Hamiltonian TL are solutions of the polynomial equation defined by Eqs.([l5|) and (|l3| ) 
through the condition det \TL — ojI = 0. To illustrate this point, figure || displays the band structure of a superlatticed 
system consisting of two materials A and B whose tight-binding parameters are arbitrarily chosen as e A = e° B = 0, 
t A = —1/2 and t B = —1/3. Two cases are shown in figure namely N = 3 and N = 6. In each case a number of 
possible values for Na and N B is presented. Due to the equivalence between the single-band-cellular and multi-band 
linear Hamiltonians, it is clear in all the cases that, for a given Na, there are N solutions per q± forming the respective 
bands. All bands are separated by energy gaps except for the ones in figure ||A corresponding to Na = and Na = 3. 
Those two cases correspond to single materials B and A, respectively. 

Since we are interested in determining the oscillation periods of the one-electron propagators, what is needed 
from the polynomial equation defined by Eqs.(|l5|) and ( |j"3| ) is the wave vector q± with which electrons with energy u> 
propagate across the superlatticed spacer, that is, the poles q± of the integrand in Eq.(|ll|). Whereas the determination 
of the eigenvalues requires finding the zeros of an iV-th degree polynomial in to, the calculation of the wave vector 
qj_ for a given energy w is straightforward due to the simple dependence of Eq. ( |l5| ) on q± . The advantage of writing 
the expression of the determinant in terms of the number of planes Na and N B becomes clear, for the oscillation 
periods can be written in terms of those parameters. Furthermore, bearing in mind that the wave vectors q± are the 
perpendicular components of the spacer Fermi surface, we can describe how the Fermi surface is deformed with the 
composition of the cells. The wave vector q± is such that 

cos(q±Nd) = Nb {det \H A {N A ) - wi| x det \H B {N B ) - oj1\ 

(tA) (t B ) 

-2t A t B det\H A (N A -l)-uji\ x det \H B {N B - 1) - wl| 

+ t\ t% det \H A (N A - 2) - lu1\ x det \H B (N B - 2) - , (16) 

where, for simplicity, the interplane hopping tA B was assumed to be the geometrical average between t A and t B . Note 
that for a given energy u (and qy) there is only one value of cos(q±N d) defined. This indicates that the off-diagonal 
propagators oscillate as a function of the number of cells N c with a single period. In this case the Fermi surface is 
deformed as a function of N A and N B but the number of periods is left unchanged. This can be understood through 
the size reduction of the Brillouin zone along the perpendicular direction q z , folding the Fermi surface back into the 
zone when it crosses the boundaries in q z = ±^j- 



B. Deformation of the Fermi surface 



Here, given two materials A and B, we shall look at the change of the Fermi surface for different thicknesses N A 
and N B . Still assuming the tight-binding parameters used in the evaluation of figure ||, the cross sections of the 
superlatticed spacer Fermi surface for q x — are shown in figures ^ and |^. Two cases are displayed, namely N = 3 
with Ep = —1.6 in figure [| and N — 6 with Ep = —1.35 in figure |], where Ep represents the Fermi level. In all 
those figures the vertical axes are along the direction perpendicular to the layers whereas the horizontal direction 
represents q y . The Fermi surfaces are drawn inside the first Brillouin zone represented by the rectangles, which are 
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reduced along the z direction by a factor N. Note that figures |^A (Na — 0) and ||D (Na — 3) correspond to the bulk 
Fermi surfaces of the pure materials B and A, respectively. The choice of different Fermi levels for N — 3 and TV = 6 
results from the analysis of figure 0, using the criterion of selecting a Fermi level which crosses the most number of 
bands. In fact, in figure |]A Ep = — 1 .6 yields solutions for all possible values of Na and Nb whereas in figure ^3 
Ep = —1.35 intersects the curves Na — 2, Na — 3 and Na = 4. 

As discussed in section ||, in the case where only one period exists in the propagators, the RKKY geometrical 
picture coincides with the more general stationary phase conditions of Eq. (]7j) for selecting the effective oscillation 
periods of the coupling. Due to the equivalence between q x and q y , it is clear from figures | and | that all wave 
vectors q± satisfying the stationary phase conditions are located at the origin of the 2-dimensional Brillouin zone, i.e., 
qx = Qy = 0. Therefore, the oscillation period of the propagator goN a is associated with the value of q± determined 
by the solution of Eq. ([l6| ) for q x = q y = 0. 

According to Eqs. (J^Jand (||), the coupling is a function of the propagators g between the magnetic planes. In the 
case considered here, 3ojv c an d 9n c % are * ne propagators involved. We recall that in this case the coupling depends 
on the product 9qn c 9 < n°o anc ^ bearing in mind that both propagators oscillate with the same periods, the coupling 
oscillates with an effective wave vector which is twice that for g$ n ■ The corresponding periods of the coupling are 
then given by A(Na, Nb) = 7r/<7_i_. The periods depend on the intra-cell composition determined by Na and Nb as 
can be seen by the change in caliper of the extremal wave vectors spanning the Fermi surface in the centre of the 
Brillouin zones. Following Eq.([l6|), the periods are displayed in Table | and the product <7jv" 9on c as a function of 
N c for different intracell parameters and for a = a' = 1 is shown in figure ||. Note the agreement between the periods 
in figure || and the ones in Table |. 



IV. DISCUSSION AND CONCLUSIONS 



The results presented above indicate that the oscillation periods of the interlayer coupling across supcrlatticed 
spacers are related to the shape of the Fermi surface when the coupling is investigated as a function of the number 
of cells of the superlattice. Therefore, a change in the composition of the superlattice enables one to deform the 
corresponding Fermi surface and consequently tune the oscillation periods of the coupling. Moreover, although the 
present analysis is not concerned with the strength of the oscillation periods, it should be mentioned that the change 
in intra-cell parameters not only alters the relevant wave vectors but also the respective curvatures of the Fermi 
surface. Bearing in mind the dependence of the oscillation amplitude on such curvatures, the contribution of periods 
may be enhanced by the variation of the parameters Na and Nb- Non-RKKY periods, for instance, which in some 
cases contribute weakly to the oscillationsEm, could become more relevant through the variation of those parameters. 

It is hoped that the simple case treated here reflects what occurs in real systems, which may be confirmed by 
future experiments. Alloys that, under certain critical temperatures, present ordered phases may be used to test the 
veracity of the present claims. In the ordered phase, Cu-Zn alloys, for instance, have a BCC structure consisting of 
interpenetrating simple cubic sublatticesEZI. In this case, cells formed by one plane of Cu and one plane of Zn are 
periodically repeated, corresponding to a possible realization of the system proposed in this communication. The 
lattice parameters of the ordered Cu-Zn alloy and that of Fe are similar, which may allow them to be grown in layered 
structures. Although the oscillation periods of the coupling cannot be tuned in this system, it is suggested that the 
periodicities associated with the corresponding Fermi surface will become evident only when plotted as a function of 
even (or odd) number of spacer ] g y prQ Alternatively, by growing superlattices of non-magnetic metals, e.g. noble 
metals, in wedge-shaped systemsEj'tEl one can change the intracell composition of the system to test whether the 
periods of the coupling can be tuned, as the present calculations indicate. 
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FIG. 1. Schematic representation of a superlatticed spacer. 



(A) (B) 




q (ji/Nd) q(7t/Nd) 

FIG. 2. Band structures E a (q±,q x — q y — 0) for (A) N = 3 and (B) iV = 6 as a function of qx- In figure A the dotted line 
is Na = 0, the full line is Na = 1, the dashed line is Na = 2 and the dot-dashed line is Na = 3. In figure B the the full line is 
Na = 2, the dashed line is Na = 3 and the dot-dashed line is Na = 4. 
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FIG. 3. Fermi surfaces for N = 3. (A) N A = 0, (B) Na = 1, (C) Na = 2 and (D) 7V A = 3. Note that figures (A) and (D) 
correspond to the bulk Fermi surface of the single-material spacers B and A, respectively. The rectangles are the Brillouin 
zones, whose dimensions are 2-zr /d horizontally and 2n/3d vertically. 



10 





FIG. 4. Fermi surfaces for N = 6. (A) Na = 2, (B) Na = 3 and (C) Na = 4. The rectangles are the Brillouin zones, whose 
dimensions are 2n/d horizontally and 27r/6d vertically. 
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FIG. 5. 

9n"o 9on* c ( m arbitrary units) as a function of the number of cells N c for some intracell parameters, a = a' = 1. 
The full lines with squares are the real part of the functions whereas the dashed lines with circles are the imaginary part. The 
symbols represent the integer values of N c , whereas the lines are a guide to the eye. (A) Na = 1 and Nb = 2; (B) Na = 2 
and N B = 1; (C) iV^ = 3 and Nb = 3; (D) N A = 4 and Nb = 2. In cases (A) and (B) the Fermi level is E F = -1.6 whereas 
E F = -1.35 in cases (C) and (D). 
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N = 3 


N = 6 


A(0,3) = 1.13 


A(2,4) = 1.69 


A(l,2) = 1.37 


A(3,3) = 3.56 


A(2, 1) = 3.67 


A(4, 2) = 8.39 


A(3,0) = 9.33 





TABLE I. Oscillation periods A(Na,N b ) (in units of D) of the coupling as a function of the number of cells N c for different 
values of Na and Nb ■ 
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